Von Koch and Thue-Morse revisited 



J. -P. Allouche 
CNRS, LRI 
Batiment 490 
F-91405 Orsay Cedex (France) 



G. Skordev 
CEVIS, Universitat Bremen 

Universitatsallee 29 
D-28359 Bremen (Germany) 
skordev@cevis . uni-bremen . de 



allouche@lri . f r 



Abstract 



We revisit the relation between the von Koch curve and the Thue-Morse sequence 
given in a recent paper of Ma and Goldener by relating their study to papers written by 
Coquet and Dekking at the beginning of the 80s. We also emphasize that more general 
links between fractal objects and automatic sequences can be found in the literature. 

1 Introduction 

The recent paper of Ma and Holdener [27] gives an interesting explicit relation between the 
von Koch curve and the Thue-Morse sequence. We show here that such a link was already 
noted in two papers dated 1982-1983, one by Coquet [8] (where the name "von Koch" is 
not explicitly mentioned, but a construction a la von Koch is given [8j p. Ill]), the other by 
Dekking, see [HI p. 32-05 and p. 32-06]. We explain the similarity of the points of view of 
these three papers through the approach of Dekking in [9] (where the Thue-Morse sequence 
does not appear). Actually the fundamental object is the sequence ((— l) S2 ^j n ) n >o and its 
summatory function, where S2{n) is the sum of the binary digits of the integer n, hence 
((-l)«W) n > is the ±1 Thue-Morse sequence, and j = e 2i7r / 3 . The sum £ n<JV (-l) aa(n) .7 B 
was explicitly introduced by Coquet in [8] to compute X/n<Jv( — -0 anc ^ m a more general 
form by Dekking in [11]. Joining the points of the complex plane that correspond to the 
sums E n<A r(— l) S2 ( n )j n gives a non-renormalized version of the von Koch curve. We end this 
paper with a quick survey of some general relations between automatic sequences and fractal 
objects. 

We suppose that the reader is familiar with the construction of the von Koch curve (see 
in particular the original paper of von Koch [26] or the book of Mandelbrot [28]) and with 
the Thue-Morse sequence (see for example [3]). We will nevertheless recall the definition of 
the Thue-Morse sequence as well as the definition of morphisms of the free monoid and of 
automatic sequences: for more details the reader can look at [I] for example. 
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2 A few definitions from combinatorics on words 



In what follows an alphabet is a nonempty finite set. Its elements are called letters. A word on 
the alphabet A is a finite string of elements of A. The length of a word is the number of letters 
composing this word. The concatenation of two words v = aia 2 ■ ■ -cij and w = b\b 2 ■ ■ -bj, 
where the a's and 6's are letters is noted vw and defined by vw := a\a 2 ■ ■ ■ Ojb^b^ ■ ■ ■ bj. The 
free monoid generated by A is the set of words on A, including the empty word, equipped 
with the concatenation of words: it is denoted by A*. 

Definition 1 Let A be an alphabet. 

A morphism ip on A is a map from A* to itself that is a homomorphism for the concate- 
nation. (The morphism ip is determined by the images (p(a) for a £ A.) 

A morphism (p on A is called d-uniform for some integer d > 1 (or uniform of length d 
or uniform if d is clear from the context) if the words (p(a) have length d for all a G A. 

Example Let A := {0, 1}. Define the 2-uniform <p on A by (p(0) := 01, <£>(1) := 10. Then, 
denoting by <p> k the fcth iterate of (p, we have: 

y?(0) = 01, <^ 2 (0) = 0110, <^ 3 (0) = 01101001, ... 

It is easy to see that the sequence of words ({p k (0)) k simply converges to the infinite sequence 

0110100110010110 ... 

called the Thue-Morse sequence on the alphabet {0, 1}, and that, for n > 0, the nth term 
of this sequence is equal to the parity of s 2 (n), the sum of the binary digits of the integer n. 
Replacing by a and 1 by b yields the Thue-Morse sequence on the alphabet {a, b} (called 
the ±1 Thue-Morse sequence if a = 1 and b = — 1) 

abbabaabbaababba ... 

Definition 2 Let A be an alphabet and d be a positive integer. An infinite sequence (u n ) n >o 
on A is called (i-automatic if there exist an alphabet B, a letter b e B, a <i-uniform morphism 
ip> on B and a map f : B ^ A such that the sequence of words (ip k {b))k simply converges to 
an infinite sequence (v n ) n >o on B, and such that for all n > we have u n = f{y n ). 

3 The result of Ma and Holdener 

The authors of [27J p. 202] describe the picture obtained with a "turtle program" (see [1] for 
example, see also[25]) based on the Thue-Morse sequence. After scaling, the picture converges 
to the von Koch snowflake. Up to designing a "simplified" curve, the authors prove in fact 
that the classical iterative construction scheme for the von Koch curve is precisely given by 
the Thue-Morse sequence in this setting. We begin with a simplified formulation of their 
result. 
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Theorem (Ma, Holdener |27| ) Let {F, L} be the turtle alphabet where F is move one 
unit forward and L is rotate by ir/3. Let TMj, be the prefix of length 2 k of the Thue- 
Morse sequence on the alphabet {F, L}. Let W n := K 2n (TM in ) be defined as follows: the 
word TM^n is written as the concatenation of words from {TM2 n ,TM2 n }, where TM.2 n is 
obtained from TM.2 n by interchanging the symbols F and L; this concatenation is interpreted 
as a word on the turtle alphabet {F, L} thus defining a polygonal line in the plane denoted 
by W n . Let S n be the scaling factor defined by S n := 2/(3 n ~ 2 — 1). Then the sequence of 
polygonal lines S2 n W n converges to the von Koch curve. 

Proof. This is essentially Theorem 5.0.14 of [27] . □ 

One sees that, up to scaling, the polygonal line W n+ i is obtained from W n in four steps: 
draw W n , glue a copy of W n rotated by — n/3, glue a copy of W n rotated by — n/3 + 2n/3 = 
tt/3, and glue a last copy of W n rotated by — n/3 + 2n/3 — 7r/3 = 0. In other words, the 
polygonal line W n is transformed into W n +i exactly as follows, up to scaling: each segment in 
W n is interpreted as a complex number z in the coordinates where the origin is an endpoint 
of this segment and the x-axis is parallel to this segment; this complex number is then 
replaced by the four complex numbers z, —jz, —j 2 z, z (where j := e 2tn ^ 3 ); the segments 
corresponding to z, —jz, —j 2 z, z in the same coordinates are glued consecutively, i.e., the 
segment we started from is replaced by the segments corresponding to the complex numbers 
z,z- jz, z- jz- j 2 z, z- jz- j 2 Z + z. 

In other words, still ignoring scaling, the curve which is obtained is the geometric repre- 
sentation of the sum J2o<k<4 n u k where the sequence (w n ) n >o is the fixed point beginning in 
1 of the morphism defined over the alphabet {1, — 1, j, —j 2 } by 



A : ( 



1 - 


* 1 


-J 


-J 2 


1 


-1 - 


* -1 


j 


3 2 


-1 


j - 


* j 


-J 2 


-1 


3 


-J - 


* -J 


J 2 


1 


-3 


J 2 - 


* J 2 


-1 


-j 


J 2 


-J 2 - 


- -J 2 


1 


J 


-J 2 



which easily implies by induction on n that 

where S2{n) is the sum of the binary digits of the integer n, and where thus the sequence 
((— l) S2 W) n > is the ±1 Thue-Morse sequence. Note that, to recover the usual (bounded) 
von Koch curve, the appropriate scaling here is 3~ n . 

This yields two consequences: 

• The occurrence of the Thue-Morse sequence in [27J parallels its occurrence in [5], where 
Coquet, interested in the behavior of the summatory function ^ fc<n (— l) S2 ^ 3k \ intro- 
duced the sum 

2~2k<n(- 1 Y 2{k) 3 k - He obtains (page 111) by-product the von Koch 
curve, which he calls "a fractal object built on a classical fractal scheme" (see also 

mi). 
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The morphism above is a simplification of the morphism given by Dekking in [9] to 
construct the von Koch curve. Namely, Dekking considered the morphism defined on 
the 8-letter alphabet {s, n, o, w, s', n', d, w'} by 
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followed by the map s — > 1, w — > 1, n — > —j 2 , o — > —j, s' — > —1, w' — >■ —1, nl — > j 2 , 
d — > j (actually there is a misprint in Dekking's paper for the images of o and d by 
this map). Dekking's morphism clearly gives the same fixed point beginning with s as 
the morphism 
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followed by the map s — > 1, n — > — j , o — > — j, s' — > —1, n' — > j , o' — > j. This last 
map being a bijection from {s,n,o,s',n',o'} to {1, — 1, j, — j, j 2 , — j 2 }, we can write 
this morphism directly on the alphabet {1, — 1, j, — j, j 2 , — j 2 } without changing the 
final sequence: 
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(this morphism is the morphism A we obtained above up to changing j into y 



4 A generalization 

Let p and q be two integers > 2. Define s p (n) to be the sum of digits in the p-ary expansion 
of the integer n. Define r] p := e 2l7r//p and r] g = e 2m ^ q . The general sum 



Z(N,p,q):= 



0<k<N 
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was considered by Dekking in [TT]. In the case where p = 2, hence r\ v = — 1 and s p = S2, the 
same procedure as in [HUH] can be applied, replacing the above morphism of length 4 

z — > z — jz — j 2 z z 

by the morphism of length Q 

Z -> (_l)«(l)^ z (_1)^(2)^ ... (-IJ-aW-l)^-^ 

where Q := 2^ 9 > and </j is the Euler function. All details are given in [TT], where Dekking 
gives in particular (see [TTJ Fig. 1, p. 32-06]) fractal pictures associated with p = 2 and 
q G {3, 5, 7, 13}: the case q = 3, i.e., Q = 4 and ?7 g = j, shows of course the von Koch curve. 

5 Automatic sequences and fractal objects 

The purpose of this section is to emphasize that there is a general link between automatic 
sequences and fractal objects. A very first remark is that being a fixed point of a (uniform) 
morphism means having some sort of self-similarity: since an automatic sequence is infinite, 
a rescaling process is often necessary to "see" the fractal object. Let us begin with a very 
simple example. Define the morphism of length 3 on the alphabet {0, 1} by 

^ 000 1 -> 101. 

Its fixed point beginning in 1 is obtained as the limit of the following sequence of words 

1 

1 1 

10 10 10 1 

101000101000000000101000101 

The reader will have immediately noted that, starting with 101, representing l's by black 
horizontal segments of length 1/3 and 0's by white horizontal segments of length 1/3, and 
"renormalizing" by a factor of 1/3 at each step yields the classical ternary Cantor set. This 
approach can be generalized to situations where the representation does not necessarily lie 
on a straight line: see in particular the papers of Golomb [22] and Giles [Till [201 I2T] where 
fixed points of morphisms and automatic sequences are not mentioned but can be found, 
and the papers of Dekking 0, QUI [TTJ . In the same direction see the papers of Dekking, 
Dekking and Mendes France, Deshouillers, Mendes France, and Mendes France and Shallit 
[T2l [T3l \15[ [29| [30] . but also the paper of Salon with an appendix by Shallit [32], where 
two-dimensional automatic sequences enter the picture, without forgetting the article by 
Prusinkiewicz and Hammel [21], and the article by Allouche, Allouche and Shallit [2]. We 
would also like to quote two papers by von Haeseler, Peitgen, Skordev [23], [21] and the 
references therein, and a paper by Barbe and von Haeseler [5]. 

Another relation between fractal objects and fixed points of morphisms or automatic 
sequences (a n ) n > occurs in the study of the summatory function ^2 n<x a n , whose asymp- 
totic behavior often involves a logarithmically periodic, continuous and nowhere (or almost 
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nowhere) differentiable function: this non-regular behavior corresponds to a geometric object 
which is fractal. A nice and usually difficult question is to determine the Holder exponents 
of the non-regular terms in the asymptotic behavior of these summatory functions. While 
the first such study, due to Delange [H], involves the summatory function of the sequence 
(s2(n)) n >o, several papers deal with the summatory function of fixed points of morphisms or 
automatic sequences, or even with more general classes of sequences: see in particular the 
work of Dumont, and Dumont and Thomas, in particular [TBI [T7] , the thesis of Cateland [6J, 
the paper by Flajolet, Grabner, Kirschenhofer, Prodinger, and Tichy [18] . and finally the 
paper of Tenenbaum [33] . 

A third direction we want to mention concerns noiselets with in particular an occurrence 
of the (2-automatic) Shapiro- Rudin sequence, see [TJ. 

What precedes should not hide the specificity of the case addressed here, namely the close 
relation between the von Koch curve and the Thue-Morse sequence. While the relations 
explained in this section are quite general, they either link a well-known fractal object to an 
ad hoc automatic sequence, or a well-known automatic sequence gives rise to some unnamed 
and not "classical" fractal object. In the case of von Koch and Thue-Morse, the fractal and 
the sequence were studied quite independently, before a link was found more then eighty 
years after von Koch's paper and more than seventy years after Thue's paper. 

Acknowledgments This paper was written during two stays at Cevis-Mevis of JPA who 
thanks heartily all the colleagues in Bremen, and in particular H.-O. Peitgen, for their very 
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